Draft version September 24, 2012 

Preprint typeset using I^T^^X style emulatcapj v. hjljW 



LAUNCHING AND QUENCHING OF BLACK HOLE RELATIVISTIC JETS 
AT LOW ACCRETION RATE 

HUNG-Yl Pu\ KOUICHI HlROTANI^, AND HSIANG-KuANG ChANG^'^ 
Draft version September 24, 2012 

ABSTRACT 

Relativistic jets are launched from black hole (BH) X-ray binaries and active galactic nuclei when the 
disk accretion rate is below a certain limit (i.e., when the ratio of the accretion rate to the Eddingtion 
accretion rate, m, is below about 0.01) but quenched when above. We propose a new paradigm to 
explain this observed coupling between the jet and the accretion disk by investigating the extraction 
of the rotational energy of a BH when it is surrounded by different types of accretion disk. At 
low accretion rates (e.g., when m < 0.1), the accretion near the event horizon is quasi-spherical. 
The accreting plasmas fall onto the event horizon in a wide range of latitudes, breaking down the 
force- free approximation near the horizon. To incorporate the plasma inertia effect, we consider the 
magnetohydrodynamical (MHD) extraction of the rotational energy from BHs by the accreting MHD 
fluid, as described by the MHD Penrose process. It is found that the energy extraction operates, and 
hence a relativistic jet is launched, preferentially when the accretion disk consists of an outer Shakura- 
Sunyaev disk (SSD) and an inner advection-dominated accretion flow. When the entire accretion disk 
type changes into an SSD, the jet is quenched because the plasmas brings more rest-mass energy than 
what is extracted from the hole electromagnetically to stop the extraction. Several other observed BH 
disk-jet couplings, such as why the radio luminosity increases with increasing X-ray luminosity until 
the radio emission drops, are also explained. 

Subject headings: accretion, accretion disks — black hole physics — Galaxies: active — magnetic 
fields — MHD — X-rays: binaries 



1. INTRODUCTION 

Black hole (BH) relativistic jets are observed from black hole X-ray binaries (BHXBs) and active galactic nu- 
clei(AGNs). A puzzling coupling between the jet and the accretion disk has been recognized. Observationally, 
the radio luminosity (which is likely due to the radiation from the relativistic jet or may also be from the semi- 
relativistic disk wind) increases with increasing X-ray luminosity (which is likely due to the radiation from the disk), 
until it shows a sudden drop when the X-ra y luminosity exceeds a certain limit (jCallo et al.l 120031 : iMaccarone et al.l 
l2003t iFender et al.ll20(Jl : iTrump et al.ll20lTI ). In terms of the dimensionless accretion rate to, this limit has a value 
about ~ 0.01, where rh = M/M^dd, M is the accretion rate, M^dd — ^Edd/xc^ is the Eddingtion accretion rate, 
^Edd = 1-38 X 1O^®(M/M0)ergs"^ is the Eddingtion luminosity, x is the efficiency converting the energy of the ac- 
creting mass into the radiation energy, and c is the speed of light. In addition, when the X-ray luminosity changes, 
the X-ray spe ctral state also ch anges , indicating that the accretion disk chang es its type with the accretion rate 
(|Fender et al.ll200l[lYmnp et al.l[201lL lEsin et al]IT997HRmrillard fc McClintockl [20061 ). Since it is believed that BH 
relativistic jets are originated from t he rotating BHs that are th readed by large-scale magnetic lines and are powered 
by the rotational energy of the BH ([Blandford fc Znaiekl 119771) . an investigation of whether the BH energy can be 
extracted by hole-threading magnetic field lines when the BH is surrounded by different type of accretion disks, may 
offers the key to an understanding of the disk-jet coupling of BHXBs and AGNs. 

How the disk type varies to is described below. When the accretion rate is extremely low, the energy transfer due 
to Coulomb collisions between ions and electrons is very inefficient. As a result, the ions cannot efficien tly transfer the 
heat to the electrons, and the disk becomes an a dvection-dominated accretion flow or ADAF for short (jNaravan fc Yil 
119951 : iNaravan et"aLlll997l : lNaravan et al. |[T998D , which corresponds to a hot, radiatively inefficient, optically thin and 
geometrically thick disk. As the accretion rate increases, the outer part o f the disk cools down due to efficient Coulomb 
collisions, and becomes a cold, optically thick and geometrically thin disk ('Shakur a fc Sunvaevll973l : lNovikov fc Thornd 
[1973), hereafter we use the term "SSD " (Shakura-Sunyaev disk) to refer to such disk; consequently, a combined disk 
is formed with an inner ADAF and an outer SSD. In a combined disk, the radius where the transit ion from thin 
disk to ADAF occurs is called the transition radius, i?tr- In general, i?tr decreases with increasing m (jHonm al fl996l : 
iManmoto et al.ll2000f ). When the accretion rate further increases and exceeds a critical value with to ~ 0.01, the 
Coulomb interaction becomes so efficient at every radius of the disk that the whole disk cools down, changing into 
the SSD. The plasma rotational motion is nearly Keplerian in a SSD, whereas it becomes sub-Keplerian in an ADAF, 
because the pressure gradient force contributed in the radial force balance in the latter type of accretion disk. 

It is interesting to note that the accretion flow geometries near the horizon are quasi- spherical for all the above 
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disks (see Section 2 for more details). As a result, there are considerable plasma flow along the hole-threading line, 
reducing the electromagnetic extraction of BH energy by contributing their rest-mass energy. This consideration 
mainly differentiates our model from previous explanations of the disk-jet couplings (e.g., Meier 2005; Ferreira et al. 
2006). We emphasize that this effect is important at low accretion rate, e.g., when m < 0.1, due to the quasi-spherica l 
geometry near the event horizon; at higher accretion rate , e.g., when rh as large as > 0.3 (Abramowi cz et al1 l2010^. 
the accretion disk type further varies to a slim disk (jAbramowicz et all 119881 120101 : iSadowski 2009) and a disk-type 
geometry near the horizon is realized. Therefore, a pure electromagnetic extraction of BH energy at most of the 
latitudes becomes a good approximation. 

To consider the contribution of the quasi-spherical accreting plas ma near the BH, we adopt the magnetohydro- 
dynamics (MHD) theory and consider the "MHD Penrose process" (iTakahashi et al.lll990f ). which includes both the 
electromagnetic and plasma inertia effects, to examine the extraction of the BH rotational energ}!^. By constraining 
the magnetic field strength and the plasma density from the theory of accretion disks in a strong gravitational field, 
we investigate the importance of plasma inertia for different accretion disk types at low accretion rates. It is found 
that the MHD extraction of BH energy takes place preferentially in a combined disk and therefore enables the launch 
of a relativistic jet; however, such relativistic jets are quenched if a magnetic dominance breaks down when the entire 
accretion disk become a thin disk. We also propose a new paradigm of the coupling between the accretion disk and 
the jet for both BHXBs and AGNs. The observed launching and quenching of relativistic jets at low accretion rates 
and several other puzzling observational results, are naturally explained. 

We introduce our model in Section 2. In Section 3, we summarize necessary formulae and concepts of the general 
relativistic ideal MHD fiow in a stationary and axisymmetric BH magnetosphere. Next, the result is given in Section 
4. Finally, discussion and conclusion are respectively presented in Section 5 and 6. 

2. THE MODEL 

To qualitatively investigate whether a relativistic jet can be launched when the BH is surrounded by different type 
of accretion disks, we consider a model with following assumptions. Comments on the assumptions are provided in 
Section 5.3. 

1. Accretion disks are described by the analytical solutions adopted from the accretion disk theories. The accretion 
disk type is important since both the large-scale field strength and the mass flux per magnetic flux tube (hereafter, 
mass loading) are related to the disk. For simplicity, we adopt the analytical solutions of the ADAF and the SSD to 
represent the disk properties. For a combined disk, these two types of disk solutions are connected at the transition 
radius. The analytical solution we adopted is provided in the Appendix B. The 'middle region' solution of a SSD, in 
which the pressure and the opacity is respectively dominated by the gas pressure and the electron-scattering opacity, 
is used. 

2. Parabolic large-scale hole-threading magnetic field lines are considered. We consider large-scale magnetic field 
lines that are dragged in from a distant region by the accretion flow. It is expected that the field lines can be finally 
accumulated near the BH with a parabolic geometry. Thus, instead of solving for the field geometry, we adopt the 
paraboloidal field solution given in .Blandford fc Znaiekl (|1977l ) for the magnetic flux function A^j,, 

A^ = y {r(l -cos^) -H2(H-cos6l) [1 - In (1 + cos 6i)]} , (1) 

where the field strength is to be determined (see Appendix C for more details). Here, we focus on magnetic field 
lines that thread the horizon, since they are responsible for the extraction of the rotational energy from a rotating BH. 

S. The accretion flow has a quasi- spherical geometry near the event horizon. An accretion fiow will have a quasi- 
spherical geo metry near the horizon if th e flow become transonic before the radius of the innermost stable circular 
orbit, riscOj (jAbramowicz fc Zurek|[l981|) . Compare to a disk-like geometry, in which the accreting plasma enters the 
BH through a plunging region near the equatorial plane, the accreting plasma in a quasi-spherical geometry can fall 
onto the horizon at a wider range of latitudes. Since th e sonic point of the accretion flow is well outside of nsco for a 
typical ADAF, which has a > 0.2 (jNaravan et al.llT997D . and the sonic point is close to risco for a-n SSD, the plasma 
distribution near the horizon should have a quasi-spherical geometry for all the disk types we considered here, ADAF, 
combined disk, and SSD. Note that, for the SSD, the plasmas are conventionally supposed to plunge into the horizon 
in a disk-like geometry, which forms a striking contrast with our current picture. 

4-. The magnetic field strength is parameterized by the gravitational binding energy. Providing that the diffuse effect 
of the large-scale magnetic fleld is relatively unimportant, large-scale magnetic flcl d lines can be 'arrested' by the 
accretion disk, being dragged fro m the outer region toward the central BH ( Naravan et al.|[2003t I Spruit fc Uzdenskvl 
120051 : iRothstein fc Lovelacell2008f ). It is convenient to parameterize its strength B at radius R by (see Appendix A) 

B = B(e,i?,CT(i?)) = £V27rGMBW^: (2) 

where G is the gravitational constant, Mbh is the mass of the central BH , a is the surface density of the disk, and 
is the ratio of the gravitational binding energy of the disk at R to the large-scale magnetic field energy inside radius 
R. The determination of e is not yet clear; however, it can be related to the ionization degree of the accretion flow 
and the magnetic Prandtl number. 



In comparison, in the magnetically dominated limit, the Blandford-Znajek process IIBlandt'ord fc Znaiekl [T977I) becomes a good 
approximation of the MHD Penrose process. 
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5. The accretion disk type is a function of the accretion rate. As described in Section 1, from low to high to, the disk 
type varies as follows: (I) ADAF, (II) combined disk that consists of an inner ADAF and outer SSD, and (HI) SSD. 
Hereafter, we call the corresponding range of to as "Range I" , "Range H" , and "Range HI" , resp ectively . Many studies 
on th e transition radius i?tr between inner ADAF and outer SSD give similar results (jHonma 1996; Manmoto et al.l 
120001). To qualitati vely describe how the disk configuration changes with to , we adopt the analytical solution of i?tr 
from lHonmal (jl996f ) 

i?tr = 2.7648a'*TO-2i?g ifTO<TOcr, (3) 

where a is the viscosity of ADAF in the a-prescription (jShakura fc: Sunvaevlll973l) . TOcr — M^^/ M-^dd is the critical 
accretion rate with a typical value < 0.01, beyond which there is no ADAF solution and hence the entire dis k becomes 
the SSD type. We recognize that the entire disk becomes thin disk if _Rtr < 10i?g is met (Honma 1996; Manr noto et al.l 
|2000( ). and define the critical accretion rate, M = Mcr, such that i?tr — 10i?g is obtained at M — A/cr- By further 
assuming that the accretion disk has a finite outer radius i?out = 1000i?g, different ranges can be defined by (see also 
Figure 1) Range I: > lOOOEg; Range II: 1000i?g > i?tr > 10i?g; and Range HI: i?tr < lOiig. 
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Fig. 1. — Disk type as a function of the accretion rate m in our model. The solid curve shows the computed value of Equation (3) 
with a = 0.1. In our calculation, we assume the accretion disk has a finite outer radius iJout = 1000-Rg and assume Rtr = 10-Rg when the 
accretion rate fn equals to the critical accretion rate rhcr. Therefore, when m has a value within Range I/II/III, the disk type is of a(n) 
ADAF/combined disk/Shakura-Sunyaev disk (SSD), respectively. See the text for more details. 

Our strategy to investigate the formation of relativistic jet at different disk types is described below. In contrast 
to that there is only out flow along the large-scale, disk-threading field lines for disk winds (Blandford fc Pavnelll982l : 
IS^dowski &: Siko"rall2010t) . there are both inflow and outflow along the hole-threading field lines. The inflow (or outflow) 
is launched when the gravitational force acting on the plasma is larger (or smaller) than the magnetocentrifugal forces 
do. Such a difference allows us to monitor the outflow properties by the inflow behavior; for an MHD flow along a hole- 
threading line, a powerful outflow is realized only when the rotational energy is extracted by the inflow. Therefore, the 
condition for jet formation is to extract the rotational energy of BH, which can be examined by solving the equation 
of motion, namely, that the relativistic Bernoulli equation (BE), of the inflow along the field line (see Section 3 for 
more details). 

The influences of the disk type lie in how they affect the large-scale fleld strength [B) and the mass loading. At 
the injection point of the inflow, Rmj, by denoting the superscript "ADAF" and "SSD" as the analytical solutions of 
ADAF and SSD, respectively, and assuming that the same value of e can be applied to both ADAF and SSD, these 
two physical properties can be estimated (see Table 1) in different region of the computation domain of to accordingly. 

Range I: ADAF type. — The large-scale field at i?inj is estimated by the ADAF surface density, that is, Ii\ ~ 
S(e,i?inj,cr^°^^) by Equation (2). The mass flux per field line at i?inj is estimated by the ADAF solution as well, 
which gives pADAF(i?i„j)uADAF(i?i„j)/Bi. 

Range II: Combined disk type. — The large-scale magnetic field at i?inj, which contains an additional term due to the 
magnetic field advected from the thin disk to the ADAF, is estimated by Bn — i?(e, i?inj, o"^'^^^) + _B(e, i?ti., cr^^'^). 

The mass flux per magnetic flux tube at i?inj is computed by p^^^^(i?inj)u^'^^^(i?inj)/-Bii, since the injection point 
is located inside the ADAF. 

Range HI: Thin disk type. — The large-scale field at i?inj is estimated by B\\\ = B[e, Rini,a^^^) and the mass flux 

per field line at i?inj is estimated by p^^°(iiinj)ii^^°(i?inj)/-Biii. 

In next section, we summarize necessary concepts of the ideal MHD flow in a BH magnetosphere, including the 
relativistic BE and the MHD Penrose process. The result of our model is provided in Section 4. 

3. MHD FLOWS AROUND A ROTATING BLACK HOLE 

Statio nary an d axisymmetric ideal MHD fl ows in a rotating BH mag netosphere have been investigated in several 
studies ([Camenzind , 1986a . b.. ,1987; .Takahashi et al.l [l990; .Hirotani et al..,1992, ). In this paper, we adopt the Kerr 
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TABLE 1 

Estimating the Field Strength and the Mass Loading at Different Accretion Ranges 



Range Disk Type Field Strength^ at the Injection Point Mass Loading'' at the Injection Point 



I 


ADAF 


Bl = Binj(e,iJinj,<7ADAF(^_)) 


pADAF„ADAF/^j 


II 


ADAF+SSD 
(Combined) 


5„=Si„j(e,iJi,j,<7ADAF(^.^^.)) 

+St,(e,Rtr,<TSSD(iJt,)) 


pADAF„ADAF/^jj 


III 


SSD 


Bin =Bini(£,^?ini,'7SS°(i?ini)) 





" See Equation (2). 

' The mass fiux per magnetic flux tube. 

metric as the background geometry with signature (-,+,+,+) and use the geometrized units such that c = G — M — 1, 
that is, the leng t h sca le is in the unit of ^GM/c^\ Most of the derivations in this section follows the calculations in 
iTakahashi et al.l ()1990l ). Note that a different signature (+,-,-,-) is used in their paper. 

3.1. The Relativistic Bernoulli Equation 
In the Boyer-Lindquist coordinate, the metric is 

,n A-a^sin^6',n 4arsin^0,,, ^sin^6',,, E,, ^ 

ds^ = dt^ dtd(t) H — d<f)^ + —dr^ + Y.d6^ , (4) 

2-i 2-i 2^ ZA 

where a = J, A = — 2r + a^, S = + cos^ 9 and A= {r^ + a^Y + ^^"^ sin^ 9. The ideal MHD condition requires 
that the electric field vanish in the fluid's rest frame, 

^ F^.u^ = F^.u^ = , (5) 

where F^i, — A^^^ — „ is the electromagnetic field tensor satisfying the Maxwell's equation, is the four velocity 
of the fluid, A^ is the electromagnetic vector potential, and the comma refers to the derivative. 

The fluid equation of motion Tf^^ = can be split into the two components: the 'relativistic BE' (the poloidal 
equation) that describes the flow along the field line, and the 'Grad-Shafronov equation' (the trans-field equation) that 
describes the force balance perpendicular to the field lines, where the semi-colon refers to the covariant derivative. 
The BE can be obtained by projecting the equation of the motion along the field line; however, it can alternatively 
derived from the definition of proper time, 

= -1 . (6) 
The energy-momentum tensor T'"^ — T^^^^ + T^^^^^^^^ consists of two terms. One is the electromagnetic term, 

ToZ = ^ (F^'^K - ^g'-F'^^F^p) , (7) 

and the other is the plasma term, 

T^Lrn.-iP + P)^'^''+9"'P, (8) 
where P is the pressure and p is the energy density. The proper number density n obeys the continuity equation 

(nun-,^ - , (9) 

and the relativistic specific enthalpy fi satisfies 

= -P + = ™p[l + ^^I^L(^)7-i] (10) 

7 - 1 71 7-1 TlinjTOp "inj 

where rrip denotes the rest mass energy of the proton, 7 is the adiabatic index, and the subscript "inj" is the quantity 
evaluated at the particle injection point. The assumption of an adiabatic fiow requires that the entropy along the field 
line be constant, which ensures that the term Pinj/(ninjmp) remains constant along the field line. 

There are four more conserved quantities along the field line: the angular velocity of the field line, fip , the particle 
flux per unit flux tube, r/, the total energy of the flow per particle, E, and the total angular momentum of the flow 
per particle, L. They can deflned by (cf. Equations(2.3)-(2.6) of Hirotani et al. 1992) 

^^ = Y--j^^ (11) 
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-P(90 Fct>r Fr6 

E = -fj,ut - ^V^F^' = -fxut + , (13) 

47777 47rry 

1 . ^re -Sri 



where —g = —det | 5 |= E^sin^^ , ft = u'^/u* denotes the fluid angular velocity, and = —y/-^F^^ = 
— {A /T,) sin 6 FrB the toroidal magnetic field seen by a distant static observer with four velocity u'(^g^ ~ (1,0,0,0). 
One should note that the sign of r] is defined by the signs of it*" and F$^. Although the magnetic field is defined by 

B^. = ^\/^e.^„^F«'5nLb ; (15) 
however, it is convenient to introduce the rescaled poloidal field 

AB 

Bl = B^BAigtt + gtc^^p)-^ = (A^,a^0,b), (16) 

9t4, ~ 9tt94>4> 

where A runs over the poloidal coordinates, r and 6. Combining Equations (12)-(14), we have 

pre V-9 _ {9t4> + ^Fg,p<f>)E + {gtt + ^Fgt<t>)L 

where 



(17) 



Ko = -{g4>4>Cll + 2gt^il.F + gtt) , (18) 

and the poloidal Alfven Mach number is 



= 47r/xn j = Anijr] j = 47r/zr/Vn . (19) 
The poloidal velocity of the plasma up is defined by 

Up = u^Ur + u^ue . (20) 

Substituting Equation (17) into Equations (13) and (14), we can express the energy of the fluid, —fiut, and the angular 
momentum of the fluid, /uw^, as 



{gtt + flFgt<^)iE-nFL) + M^E 



{gt^ + nFg^^){E-nFL)-M^L 
^"^ = M2^^ • (22) 

Finally, combining Equations (4), (20), (21), and (22), we obtain the BE (cf. Equation (17) of Takahashi et al. 1990) , 

2 , , _,E^^ KoK2-2K2M^-K4M^ 
+ ^ ~ W (M2 - Kor ' ^ ' 



where K2 and K4 are deflned as 



K^^il-n^-] , (24) 



94>4> + '^9t4>^ + gtt\^ 



/{gls - 9ttg>j»j>) • (25) 



The BE, Equation (23), can be rewritten as a polynomial of up with the aid of Equations (10) and (19). Note that the 
order of the polynomial equation depends on the second term of Equation (10), because particle number conservation 
gives n oc 1/up. 
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3.2. The Cold Limit 

The cold hmit can be adopted when the second term in Equation (10) is relatively unimportant to its first term, 
i.e., -Pinj/'T'inj^pC^ (in c.g.s unit). Equation (10) therefore reduces to ^ = and the BE becomes a fourth-order 
polynomial equation of up. By introducing a parameter C, = Up/Bp, we have 

Asm 
^ Asm 

The cold BE, therefore, can be expressed in terms of ( 

{A^B)C^ - {2KoAB)C^ + {A^ - A^C + K^B)C^ + {2VA ~ 2KaA)C + K^^K^ -V)=Q, (26) 

where 

A = 47r/x77 , 
B= =Bl , 

M 

Even for a disk whose temperature is comparable with the virial temperature and the internal energy is comparable 
to the released gravitational energy (e.g., ADAF), the cold limit can still be modestly satisfied. This is because the 
released gravitational binding energy is still a small fraction of the rest-mass energy. Therefore, we solve the BE in 
cold limit. Equation (26), for simplicity and expect the result is at least qualitatively correct. In order to specify the 
coefficients A and S, we need to specify the mass loading (mass flux per magnetic flux tube), /xry, and the strength 
and geometry of the magnetic field. It is provided in the Appendix C how we compute these coefficients. 

3.3. Light Surfaces and the Separation Point 

There are two light surfaces, which are defined by Kq = 0, in a BH magnetosphere. The outer light surface is formed 
by the centrifugal force in the same manner as the light cylinder in a pulsar magnetosphere, and the inner light surface 
is, on the other hand, formed by the strong gravity of the BH and the rotation of the magnetic field. 

If a plasma starts with a negligible poloidal velocity, the inflow and outflow along the field line separate at the 
point Ts where the gravitational force balances with the magnetocentrifugal forces. This point is called the "separation 

point" (or the "stagnation point") and defined by Kq = 0, which leads to (In up) — 0, where the prime denotes the 
derivative along the flow line. That is, if a plasma is injected at r = ri„j with ri^j < rg (or with rinj > Vg), it will be 
accelerated inward (or outward), subsequently passing through the Alfven point r — r^ and the inner (or outer) light 
surface. Figure 3 shows the contours of Kq around a BH with spin parameter a = 0.6. Note that, at the injection 
point, the BE reduces to 

E-nFL = fi^olr,^^ . (27) 
3.4. Critical Points 

The critical points for an MHD flow can be foun d by differentiating E qua tion (23) along the flow line, (In wp) = N/ D, 
where explicit expressions of N and D are given in lCamenzindl (|T986b) and lTakahashi et al.l (jl99 0^. For a general MHD 
inflow, there are three critical points at which both D and TV vanish: the 'slow-magnetosonic point', the 'Alfven point', 
and the 'fast-magnetosonic point'; the poloidal velocity Up there matches the slow-magnetosonic velocity, the Alfven 
velocity, and the fast-magnetosonic velocity, respectively. In the cold limit, the slow-magnetosonic velocity reduces to 
zero. Due to the causality requirement at the event horizon, any MHD inflow must pass through the fast-magnetosonic 
point, after passing through the Alfven point. 

At the Alfven point r = ta, the Mach number equals Kq, 
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4 6 

Rxsine / (GM/c^) 

Fig. 2. — Illustration of the separation point of a field line and contours of Kq around a black hole with spin parameter a = 0.6 and 
= 0.5Qh- Light surfaces are defined by Kq = 0. The Alfven point of a physical MHD flow must be located inside the two (inner and 
outer) light surfaces, because M^\rj^ = Kg |r^> 0. Black region represents the BH. For a field line threading the BH (black solid line), the 
separation point is determined by Kg = (indicated by a cross). See the text for more details. 



(28) 



Note that the requirement of Af^ 



Kq |rA^ implies that ta must be located between these two light surfaces. 



Besides, it is convenient to express E and L in terms of the quantities evaluated at the Alfven point. Since both D 
and N automatically vanish at the Alfven point, no additional constraints are imposed on the MHD flow at the Alfven 
point. Therefore, by requiring that the numerators in Equations (21)-(23) vanish at r — ta, we obtain 



E 



{gtt + ^Y9t^){E - VLpL) 
M2 

(gt0 + VL-Fg4,4,)(E - npL) 
M2 



and 



(29) 



(30) 



K2 |rA= ~K2Ko W ■ (31) 

Combining Equations (27), (29), and (30), we can express the conserved quantities E and L in terms of rinj and r^- 

After the MHD inflow pass through the Alfven point, D vanishes again in (Inwp) — N/D at the fast-magnetosonic 
point; thus, N should also vanish there. We should notice here that both D = and N = can be satisfied at the 
fast-magne tosonic point only f or a specific combination of the conserved quantities, fip, E , L, and rj. In previous 
works, e.g.. iCamenzindl (|1986aD and iTakahashi et all (|1990D . the authors gave rip, and L (or equivalently, fip, nnj, 
and ta) and search rj to select the trajectory that pass through the fast-magnetosonic point. In this work, we instead 
give rip, Hnj, and r], considering plasma injection from the innermost region of the accretion disk, and solve for ta (or 
equivalently, E). As explained next, whether the rotational energy of the BH is extracted outward can be determined 
by the location of r^. 

3.5. Negative Energy Flow 

The sign of E can be solely determined by ta via Equation (29) because it follows form Equation (27) that E—ftpL > 
holds for any MHD flow starting from a position between the two light surfaces with a vanishing poloidal velocity. 
Note that E becomes negative if {gu + ^Fgtip)\rA > . Thu s, near the BH, a spati al region called the "negative energy 
region", which satisfies {gu + i^p^t^) > 0, can he defined (jTakahashi et al.lll990[ ). If the Alfven point resides in this 
region, the MHD inflow has E < (which also implies L < because oi E — ilpL > 0; whereas L < does not imply 
E > 0). 

A spinning BH differs from a static BH by an additional structure called "ergosphere" (which is defined by gtt > 0), 
a region inside which any particle must orbit in the same direction of the hole spin. Being locating outside the horizon, 
the ergosphere can be viewed as the region that stores the rotation energy of a BH. Note that the negative energy 
region exists only when fin > ^^f > 0, and it is always located inside the ergosphere (see also Figure 3). It is also 
noteworthy that a negative energy MHD infiow can be launched from a region outside the ergosphere. 
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Fig. 3. — Spatial position of the Alfven point that gives a negative energy flow for the case of Qp = 0.5Qh and a = 0.6. The dash-dot 
line, which is defined by gtt + ^Fgt4, = 0, is located between the inner light surface (thin solid line) and the static limit (dotted line), the 
latter which defines the outer boundary of the ergosphere. An MHD flow has a negative energy if the Alfven point is located within the 
negative-energy region, which is defined as the region between the das h-dot and the thi n solid lines. The magnetic field threading the hole 
surface in Figure 2 is also shown (thick sold line). Sec also Fig ure 6 of lTakahashi et al.l ||T990|) . 



TABLE 2 

Inflow Energy as the Criterion for Ergospheric Jet Launch 



Inflow Property 




Infow Energy'' 


Inflow Energy Flux'^ 


Jet Launch?'^ 


plasma dominated {£1^ < 


plasma 1 ' 


E >0 


f < 


NO 


magnetically dominated (f^^j 


> If \) 
1 plasma 1 ' 


E <0 


> 


YES 



° fom > 3.nd fpjjjgjjjjj < arc usually satisfied. See the text for more details. 

^ The total energy E, Equation (13), contains both the electromagnetic and the plasma contribution. 
S'^ = nEu^ = + Spi^g^^, Equation (32), where u'" < holds for an inflow. 

A relativistic jet, which has > 0, E > 0, and hence f > 0, can be launched only when it connects to an inflow with £'~ > 



3.6. Outward Energy Flux and the MHD Penrose Process 
The outward energy flux of the flow f can be split into the electromagnetic part and the fluid part ^^piasma' 



'-' — — "-^ " ^ "--cm 1" <^plasma ) 



where, 



cr ^ 
'-'cm — 



ftp 

47Esin6'' 



cr 
plasma 



-nilUtU 



(32) 

(33) 
(34) 



Therefore, a negative energy inflow (E < and < 0) results in a positive outward energy flux (f" > 0), that is, the 
BH rotational energy is extracted. Note that £1^ > usually h olds because its sufficient condition, < < Qh, 
is satisfied under normal conditions ()MacDonald fc Thornel[T982D . where fin refer to the rotational angular velocities 
of the BH. Note also that the plasmas usually cannot extract BH's energy (i.e., Spiasma usually holds) because of 
their rest-mass contribution, unless they are in a negative energy orbit, which is rare. 

In a magnetically dominated magnetosphere ( > fpiasma)' because f is continuous across the separation point, 
a stationary outflow solution (u'" > 0, £' > and £^ > 0) is possible only when it connects a negative-energy inflow 
{ E < and u"^ < 0, such that f > 0). Consequently, a stationary MHD jet from a BH ergosphere ( a "ergospheric 
jet") can be launched in a magnetically dominated magnetosphere only when the MHD Penrose process, f" > , 
operates. If a magnetic dominance breaks down, can be discontinuous across the separation point, that is, an 
outflow (f > 0) can connect to a positive-energy inflow < 0), owing to the energy supplied by the accreting 
plasma near the separation point. However, this kind of outflow, similar to the disk wind, will not be accelerated into 
relativistic energies and therefore is not of interest in the present paper. To summarize, for an MHD flow along the 
hole-threading field line, a stationary outflow with positive outward energy flux (f > 0) should connect an inflow 
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Fig. 4. — Extraction of BH rotational energy by MHD inflow (left column: for a BH mass of IOMq; right column: 10^ Mq). The domain 
of accretion rate m is divided into three Ranges (I, II and III) according to the disk types. The parameter e (< 1) is the squared root of 
the ratio between the large-scale magnetic field energy and the disk gravitational binding energy; that is, e denotes the ability for a disk to 
confine the large-scale field, (a) and (e) Large-scale magnetic field strength B (solid line) at the injection point. In Ranges I and III, only 
(dashed line) contributes to B, where B^^j denotes the field strength evaluated at the injection point. In Range II, an additional term 
Bti (dash-dotted line) is added so that the fields advecting from the outer SSD into the inner ADAF region may be taken into account, 
where Btr denotes the field strength evaluated at the transition radius of the SSD and the ADAF (see Section 2 for more details). For 
comparison, we also plot the strength of local magnetic field when the plasma beta becomes 1 at the injection point as the dash-dot-dot 
line. To avoid complexity, only the case of e = 0.1 is shown here, while e = 1, 0.2, 0.1 and 0.05 arc shown in the rest of panels, (b) and 
(f) The location of the Alfven point R\. (c) and (g) Zoom-in plot of the second panel, (b) and (f), near the ergosphere. The static limit 
indicates the outer boundary of the ergosphere. The flow energy becomes negative if R\ is located inside the 'negative energy region' 
indicated by the shaded region (see also Figure 3). (d) and (h) Total energy (consisting of plasma and magnetic energies) per particle in 
unit of particle rest mass energy, rUpC^ . Unless E fa rripC^ , the magnetosphere becomes magnetically dominated. 



solution with f > at the separation point, provided that the injected energy flux by accretion is small compared to 
the Poynting flux. Thus, we can investigate the activity of a Poynting- flux-dominated jet (or an outflow) by examining 
the inflow energy, E (see also Table 2). 

It should be noted here that the term "MHD Penrose process" refers to the case of f > (rememb er E is the 'total' 
energy, consisting both the electromagnetic plus plasma energy of the flow) , as initially proposed by iTakahashi et al] 
([1990) : however, it was used to indicate f pias ma > (i.e., a negative energy orbit of the particle) in some previous 
studies (|Hirotani et al.lll992l:lKoide et al.|[200^ iSemenov et al.ll2004l : lKomissarovll2005D . The latter case is possible, but 
can only be transient (" Koide et al.l 20021: iKomissaro vll2005D . In comparison, in the original idea of the 'Penrose process' 
(jPcnrosc 196_9), only the particle contribution (f = fpiasma > 0) was considered; whereas in the Blandford-Znajek 
process (jBlandford fc Znaieldll977f ). only the electromagnetic contribution (£'' = E^^^ > 0) was considered. 



4. RESULT 
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We choose representative parameters to demonstrate the jet behavior at different disk type quahtatively. Because 
the extraction of the rotational energy of BH is mainly related to whether the BH magnetosphere is magnetically 
dominated or not, the change of the parameters will result in only minor modification of the result in our model. 

We adopt the dimensionless BH spin parameter to be 0.6 and the angular velocity of the field, ilp, to be half of the 
angular velocity of the hole, iln (i.e., fip = O.Siln = 0.0833). The MHD inflow along the field line that intersects the 
event horizon at a modest latitude (^h ^ 1.1868rad, i.e., about 68° from the pole; see Figures 2 and 3 for this field 
line) is considered, and the infiow is assumed to be launched at the separation point (that is, the injection point is 
located at the separation point). In addition, to investigate the importance of the BH mass, we consider separately 
the BHXB and AGN cases. For BHXB, we assume Mbh = IOMq, whereas for AGN, A/bh 10*^ Mq. Furthermore, 
to explore the dependence of the ability of the disk in storing the large-scale magnetic field, we adopt different values 
of e (e = 1,0.2,0.1, and 0.05). 

After starting from the injection point, the accreting MHD plasma is accelerated inward, passing thorough the 
Alfven point, Ra, and the fast-magnetosonic point, and finally enter the BH. The inflow energy E can be computed 
by imposing the fast-magnetosonic condition once the injection point, the angular velocity of the field and the mass 
loading are specifled. However, to discuss the conditions for jet launching (f > 0) and quenching (f < 0), it is 
adequate to examine the sign of E, which is uniquely specifled by the location of -Ra- Speciflcally, if i?A resides in 
the negative energy region, the MHD inflow has a negative energy {E < 0) and the MHD Penrose process operates. 
Therefore, by solving Ra from the BE for an MHD inflow onto a BH, we can investigate whether the energy extraction 
takes place as a function of m. Detailed computation steps are provided in the Appendix. 

Let us flrst consider the BHXB case with the BH mass, Mbh = lOAf0 (Figures 4(a)-(d)). Figure 4(a) shows how 
the large-scale fleld strength at the injection point varies with the accretion rate. It is important to note that in a 
combined disk (Region II), the surface density of the disk does not monotonically increase inward but has a local 
maximum at the transition radius, Rtr- Since an SSD can conflne more large-scale fleld than an ADAF because of 
its greater surface density for a flxcd e, the large-scale fleld arrested in the inner ADAF for a combined disk (solid 
line) is relatively stronger than what the ADAF alone would bring into the BH (dashed line). The locations of Ra 
for four discrete values of e are shown in Figures 4(b) and 4(c); the solid, dashed, dash-dotted, and dash-dot-dot lines 
correspond to e = 1, 0.2, 0.1, and 0.05, respectively. It follows from the zoom-in flgure (Figure 4(c)) that Ra enters 
the negative energy region (indicated by the shaded region) in Range II for e = 1, 0.2, and 0.1, indicating that the 
MHD Penrose process turns on in a combined disk. This is because, near the horizon, the relatively strong magnetic 
fleld (which is mainly provided by the outer SSD) and the relatively small mass loading (which is provided by the inner 
ADAF) result in a magnetically dominated environment. In comparison, a magnetically dominated environment near 
the horizon is relatively harder to achieved when the disk is an ADAF (Range I) or an SSD (Range HI), because it is 
the same type of disk responsible for both the conflnement of the large-scale fleld and the loading of the plasma. Such 
disk-jet coupling provides the reason why observed relativistic jet is launched when rh is near (but smaller) the value 
0.01, but quenched when rh exceeds ~ 0.01. In addition, since B increases with increasing to, E increases negatively 
(i.e., the rotational energy of BH is more efficiently extracted outwards) with increasing to once Ra enters the negative 
energy region in Range II (Figure 4(d)). This is because the BH magnetosphere becomes more magnetically dominated 
when the transition radius decreases with m, resulting from that the difference between the surface density of the outer 
SSD and inner ADAF becomes more significant. 

Next, let us consider the AGN case with Mbh — W^Mq (Figures 4(e)-(h)). Although the solutions look qualitatively 
similar, there are three major differences from the BHXB case. First, comparing Figures 4(b), (c), (f) and (g), we find 
that for an AGN, E < is realized in a wider range of to for a fixed e or in a wider range of e for a fixed m. Second, 
comparing Figures 4(d) and (h), we find that an AGN can extract more BH rotational energy per plasma particle 
than a BHXB can. Third, in an AGN, the MHD Penrose process works also in an SSD (i.e.. Range HI) for a relatively 
greater e (e.g., s > 0.2). However, the jet power decreases with increasing to in an SSD (i.e.. Range HI), in contrast 
to the case of a combined disk (i.e.. Range II). The above comparison can help to understand the observed properties 
of accreting BH in different masses (see below). 

5. DISCUSSION 

5.1. Universal Paradigm of Black Hole Disk-Jet Coupling at Low Accretion Rate 

The results in Section 4 can offer the key to an understanding of the observed connections between the disks and 
the jets for BHXBs and AGNs. A universal paradigm of disk-jet couplings on the m — E plane (c.f. Figures 4d & 
4h) is diagrammatically presented in Figure 5. Transitions of disk-jet states are indicated by the arrows for three 
representative values of e. Since the terminal speed of an MHD outflow is roughly determined by the magnetic 
dominance at the launching point, the more magnetically dominated the disk innermost region becomes (i.e., the more 
a system moves downward on the m — E plane), the greater speed an MHD outflow achieves. This discussion is valid 
for both the disk wind near the BH and the ergospheric jet. The two downward arrows in the right edge of this flgure 
indicate how the terminal speed changes for a disk wind and a relativistic jet. Note that the ergospheric jet can only 
be launched when E < (i.e., in the lower half of the diagram). 

The transitions of disk-jet sta tes for BHXB s , whic h are indicated by the solid arrows, explain the observed features 
of the BHXB disk-jet couplings. iFender et al.l ()2004D show that two type of outflows associated with the X-ray spectral 
types are required to explain the systematic properties of the radio emission of BHXBs. The flrst type is a steady 
outflow observed when the BHXB is in a hard state (HS), while the second type is a transient, ballistic, powerful 
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Fig. 5. — Universal paradigm of the disk-jet coupling for both BHXBs and AGNs. As rn increases, the disk-jet system transits along the 
solid arrows for BHXBs and along the dashed arrows for AGNs. The tracks of transit are depicted for three representative values of e, 
according to the result in Figure 4. Puzzling observational disk-jet features, such as why radio emission is quenched beyond a certain m 
and why the radio luminosity decreases (or increases) with increasing X-ray luminosity when the system is luminous (or less luminous), 
are explained. See the text for more details. 

and relativistic jet when the source is in a steep power-law (SPL) state and being switching from the HS to the 
soft state (SS). We interpret that the steady outflow (first type) is associated with the disk wind and the transient 
relativistic jet (second type) is associated with the ergospheric jet ejected for a modest or small e. Note that, as 
shown in the paradigm, when the entire disk become an SSD, not only that the ergospheric jet is turned off but also 
that the disk wind weakens, which explains the sudden drop of observed radio luminosity when the state transits to 
the SS. In addition, the energy of the jet per particle, \E\, as well as the ergospheric jet power, iErgojct = / £^dS, 
increase with increasing m in a combined disk, where the surface integral J dS is carried out on a closed surface 
(e.g., the horizon); this likely corresponds to the observational fact that both the Lorentz factor of a jet and its 
radio luminosity rapidly increase with increasing m before the jet is quenched. When the BH is surrounded by a 
combined disk and the MHD Penrose process takes place, the relati vely strong field for t he inner ADAF can result 
in a configura-tion which has been studied in previous simula tions (jNaravan et al.l 120031 : llgumenshchev et aTl 120031 : 
llgumenshchevi 120081 IPunslv et al."2009': 'McKinne v et al.ll2012| ). here we use the name magnetcially choked accretion 
flow (MCAFJ^ according to McKinney et al. (201^. Accompanying by the jet formation, the formation of an MCAF 
can also explain some observation features when the source is in SPL. For an MCAF, because the slowed down fiow can 
earn more time to release their heat energy via radiation, the efficiency for the flow to convert mass to radiation energy 
increases (Naravan ct al. 2003). This can be the reason why the disk luminosity in SPL state is higher than that in 
HS or SS (|Fender et al.,,2009 ). In addition, the high-frequency quas i-periodic oscillation (Q PO) feature observed at 
SPL state may be related to the jet-disk QPO mechanisms found bv iMcKinnev et all ()2012D . 

For AGNs, the transitions are indicated by the dashed arrows, which imply the following disk-jet couplings. First, 
an efficient MHD Penrose process (i.e., a greater \E\) is viable at a relatively smaller m in an AGN compared to a 
BHXB. Since the SSD exists only in the outer region w hen ifi is small, it is con sistent with the observational fact that 
radio-loud AGNs favor a spectrum without a UV bump (|Kawakatu et al.ir2009t ) , which possibly arises in the inner edge 
of an SSD. Second, the MHD Penrose process works more efficiently with increasing BH mass. This conclusion explains 
why AGN jets show higher Lorentz factors than BHXB jets. Moreover, it is also consistent with the observational 
fact (|Liu et al.l 2006; Chiabcrgc & Marconi 2011) that most radio- loud AGN host BHs with larger masses (~ 10^ Mq), 

^ An MCAF is formed when the accumulated magnetic field near the BH reaches the equipartition level. The strong magnetic field 
is able to support (or disrupt) the disk at certain radius Rm, where e = 1. Although slowed down by the magnetospheric barrier near 
Rjn, accret ing plasma can still accrete inward in the form of highly non axisy mmetry, irregular streams due to magnetic Rayleigh- Taylor 
instability l|Igumenshchev et al.l2003 : IPunslv et al.l2009l : IMcKinnev et al.l2012l ). We stress that the magnetically dominated magnetosphere 
developed when the BH surrounded by a combined disk can result in an MCAF or other conceptually similar accretion flow. 
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while radio-quiet AGNs host BHs with no obvious mass constraints. Third, ahhough a BHXB jet disappears when 
the disk transits from a combined disk to a SSD, an AGN jet could still be launched from the ergosphere, because the 
SSD in an AGN tends to be more magnetically dominated than in a BHXB owing to a smaller mass loading, pu^ /B, 
around a more massive BH. We interpret that these ergospheric jets ejected from the inner edge of the thin disks in 
AGNs result in (at least part of) the radio-loud quasi-stellar objects; this forms a striking contrast to BHXBs from 
which no strong or steady radio emission has been detected in their SS. Estimated from the dynamical timescale, the 
characteristic time scale for an accreting BH system to sustain the relativistic jet, t, is scaled by the central BH mass 
with the relation t (x Mbh- Therefore, with t equals to days for the relativistic jet from a BHXB with A/bh ^ IOMq, 
relativistic jets from an AGN with Mbh ^ IO^Mq may last for ^ 10^ years. 

Note also that the relation between the accretion rate, m, and the ergospheric jet power, iErgojct, inverses in Ranges 
H and HI. That is, iErgojct increases (or decreases) with increasing m, when m is relatively low (or high) and the disk 
is a combined disk (or an SSD). Provided that the X-ray luminosity indicates the disk accretion rate, these conclusions 
naturally explain the empirical relation for both BHXBs and AGNs: the radio luminosity increases with increasing X- 
ray luminosity wh en the bolometric luminosity, L , is small, i.e., when L/Lem — M / M^dd ~ 1% ([Merloni et al.ll2003l : 
iFalcke et al.ll2003 ) , while the radio luminosity decreases with increasing X-ray luminosity when some BH systems are 
luminous, i.e., when L/L-^dd ~ 10%(King et al. 2011) . 



5.2. The Estimated Speed of Ergospheric Jet 

At large distances, the spacetime is described by the Minkowski metric; thus, the Lorentz factor of the jet, F, is 
defined by F = u*. Because E is conserved quantitatively along the field lines and because the second term of Equation 
(13) becomes negligibly small compared to the first term at large distances, we obtain 

F = ^ , (35) 

where iJout denotes the total energy of the MHD outflow. 
We can further estimate the terminal jet speed, by linking i?out with the total energy of the MHD inflow, i^in, as: 

E,^, = ^^^E,^ = -\J^^^\E,,, , (36) 

J^outMout ^outMout 

where the subscript "in" and "out", respectively denotes the quantities for the inflow and the outflow. The above 
equation is obtained by considering the conservation of the outward energy flux near the separation point, provided 
that the magnetic field configuration does not change significantly there. In general, the ratio of the inflow and outflow 
particle flux, \ni-au[^/noMtu'^Q^^\, is greater than unity because only a small portion of the plasmas will escape as an 
outflow. As a result. Equations (35) and (36) give 

F > . (37) 

M 

Hence, the inverse values of the total energy shown in Figures 4(d) and (h), as shown in Figure 6, give an estimation 
of the lower limit of the ergospheric jet Forentz factor for BHXBs and AGNs. The BHXB and AGN relativistic jet 
speeds inferred from the observations, which has the typical value F > 2 for BHXBs (e.g.. Table 1 of Fender et al. 
2004) and F ^ 10 for AGNs, can be consistently explained by our model. 
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Fig. 6. — Lower limits of the Lorentz factor of ergospheric jets, F, modeled for BHXBs (a) and AGNs (b), as a function of the accretion 
rate, m. See text for more details. 



LAUNCHING AND QUENCHING OF BLACK HOLE RELATIVISTIC JET 



13 



5.3. Validity of the Model 

To investigate the reason of the launching and quenching of BH relativistic jets when m varies, we solve for stationary 
negative energy inflow solutions (which corresponds to the formation of relativistic jets) by adopting analytical solution 
of the ADAF and the SSD, assuming a spherical accretion geometry near the horizon, and estimating the field strength 
of the large-scale hole-threading magnetic field by the parameter e. The major drawback of our approach is that we 
ignore the possible disk instability when the SSD becomes radiation pressure dominated. In addition, our axial- 
symmetry and stationary solution miss the temporal and the non-axisymmetry properties of the jet formation. Also, 
the disk dynamo effect and the magnetic reconnection process can further modify the large-scale fields, although their 
effects can be in general included in the parameter e. 

For our purpose, we ignore several aspects which we think are minor factors in determining whether the jet is on or 
ojf for simplicity. Several discussions on the assumptions and the robustness of the result are provided in the following. 

We ignore the Blandford-Znajek extraction in the force-free region near the pole of the horizon in a quasi-spherical 
accretion fiow. However, such extraction is expected to be a relatively minor factor to determine the jet launching. 
The reason is that most of the hole-threading field lines have considerable plasma loading on them, and that the power 
for magnetic fields extracting the energy from a rotating BH is proportional to sin e (IBlandford fc Znaiekl[l97l . 

The most important feature of a combined disk is that the significant surface density difference of the SSD and 
the ADAF at the transition radius. Although here we describe the combined disk solution by patching the analytical 
solutions of the SSD and the ADAF at the transition radius, the above feature is expected to be true even if a more 
detailed disk solution is considered. 

There are uncertainties in determining the large-scale magnetic fields. For example, the parameter e can vary with 
the radius r, and the total flux of the large-scale disk-threading magnetic fields is not unclear. However, the field lines 
that threads the disk are expected to be redistributed when the accretion disk types vary. The disk-jet coupling features 
described by our model can be qualitatively preserved as long as the diffuse effect is relatively unimportant than the 
advection effect for the large-scale disk-threading field lines. Especially, a magnetically dominated magnetosphere near 
the horizon is still most likely realized when the BH is surrounded by a combined disk, provided that the large-scale 
magnetic field advected inward from the outer SSD is large enough for the inner ADAF. Such requirement is not 
difficult to be satisfied as the transition radius gradually decreases to a small enough radius (Figures 4(a) and (e)). 

The result presented in Figure 4 shows that a negative energy flow may form when the BH is surrounded by a 
combined disk (Range II). The conclusion can actually still be valid when more detailed pictures are considered. 
Below we discuss two possible concerns. First, we consider whether the result would change when numerical ADAF 
solutions, e.g., Narayan et al. (1997) and Gammie & Popham (1998), instead of the self-similar solutions. Equations 
(B1)-(B3), are used. It is true that the self-similar analytical solutions are not a good approximation near the hole. For 
example, the causality requires a numerical solution of u^^^^ be greater than Equation (Bl) near the BH, while the 
surface density a^^^^ be less than the self-similar solution. Equation (B3). Nevertheless, the numerator in Equation 
(C2) will not change because of the mass conservation, au Hpu oc M, provided that the disk height H remains 
unchanged. Since the magnetic field is predominantly supported by the outer thin disk, Bu is essentially given by 
B{Rtr)- As a result, neither /it; nor Bu contains significant errors if we adopt the self-similar solution of ADAF. Thus, 
the discussion given in this paper is not vulnerable for details of the ADAF solutions. Second, we discuss whether the 
estimation of the mass loading is fair when Bu = B{Rti ) + B{Rinj) ^ i3(i?inj). In this case, an MCAF is formed. The 
highly non-axisymmetry, irregular MCAF will try to avoid the str ong field ranges by the magnetic Rayleigh- Taylor 
instability and result in strong-field, low-density 'mag netic islands' (jPunslv et al.ll2009( ). As a resuh, both ^adaf ^^^^ 
yADAF (^ggy-gase in a magnetic island. Thus, even when > i?i„j, the mass loading from the magnetic island located 
at i?i„j is lower than the estimate with the self-similar solutions. Equations (Bl) and (B2). This means that the MHD 
inflow launched from the magnetic island should be even more magnetically dominated and that the energy of the 
inflow is still negative. 

6. CONCLUSIONS 

In this paper we focus on the fundamental question: Why relativistic jets from accreting BHs can be selectively 
launched under a certain accretion condition? Here we solve this problem by incorporating the plasma inertia effect, 
which have been neglected in previous studies, in conjunction with the electromagnetic effect on the extraction of BH 
energy. At low accretion rate (e.g., m < 0.1), the accretion flow has a quasi-spherical geometry near the BH event 
horizon. As a result, whether the rotational energy of the BH can be effectively extracted outward is determined by 
the inward rest-mass energy of the accreting plasmas. When the outward electromagnetic energy flux dominates, a 
relativistic, ergospheric jet can be therefore launched; on the contrary, when the inward rest-mass energy dominates, 
there is no powerful relativistic jet. By examining the MHD Penrose process for BHs surrounded by different types of 
accretion disk at low accretion rates, we conclude that the MHD inflow becomes magnetically dominated preferentially 
in a combined disk, which consists of an outer SSD and an inner ADAF, to enable the launch of an ergospheric jet. 
Such relativistic jets are quenched if a magnetic dominance breaks down when the entire disk becomes an SSD. We 
also provide a universal paradigm of the coupling between the BH accretion disks and their relativistic jets at low 
accretion rates. Several observed disk-jet connections are naturally explained. Future study on the relation between the 
accretion geometry near the event horizon and the extraction of the BH energy will be important to the understanding 
of the BH jet formation mechanisms. 
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APPENDIX 

PARAMETERIZED THE LARGE-SCALE MAGNETIC FIELD STRENGTH 

The large-scale field B stored in a disk can be estimated from the point of view of energies. The large-scale field 
energy inside the volume enclosed by a radius R should not exceeds the gravitational binding energy Eq of the disk 
near R, otherwise the field would escape by buoyancy. By considering Ep ^ AttR^B^ /2A and | Eq |~ 27rGAfBHO'i?, we 
define a ratio (< 1) between this two quantities to obtain 



B - ev/l2GMBH(T/i?2 , 

where G is the gravitational constant, Mbh is the mass of the central BH, and a is the surface density of the disk. 

Alternatively, we can estimate the large-scale fie ld B stored in a disk by a force balance. If a disk accretion is 
supported by a large-scale field at i?, as described in lNaravan et al.l (|2003[ ). the large-scale field B can be estimated by 
the force balance GA/bhO'/^^ BrBg/2'K to give B ^ y/2TrGMBHa/R? where Bg ^ Br ^ B is assumed. Similarly, a 
parameter £^(< 1) can be defined by the ratio between the gravity force and the field force; we thus obtain 

B - e^/2^^GMBli<y/R^ . 

Since the two estimations above are consistent in order of magnitude, we simply parameterize B in terms of e, Mbh, 
and (j/R^, such that 



With the scaled units, 



B = B{e, R, a{R)) = ex/2TTGMB}icr/R^ . (Al) 
ANALYTICAL SOLUTION OF THE ACCRETION DISKS 

A/bh = mMQ , 



R = rRg i?g = ^^^^ = (2.95 x lO^m) cm , 

M = mAhdd Mem = iEdd/O.lc^ = (1.39 x lO^^m) g/s = (2.2 x IQ-^m,) Mq/jt , 

we adopt the following radial velocity u, plasma density p, and surface density a solutions of an ADAF and a thin 
disk. Note that what we need is pu (because ^r] = pu/ti) to solve the BE, and that pu can be more or less correctly 
evaluated even in the innermost region by an ADAF self-similar solution without invo king a numerical solution. We 
therefore adopt ADAF self-similar solution (jNaravan fc Yi|[T995t iNaravan et al. iri998[) and obtain 

^ADAF ^ ^ X ^QiO(^ADAF)--i ^^^-1 ^ (g^) 
^ADAF ^ ^ ^ ^Q„4(^ADAF)-l^-l^--f ^^^-Z ^ (33) 
^ADAF ^ M/2^i?(-uADAF) ^ gg(^ADAF)-l^--i ^^^-2 ^ (33) 

while in the thin disk solutions (jShakura fc Sunvaevlll973l : iKato et al.ll2008[ ). we obtain 

^SSD ^ 20 X {a'^'^°)-Tom'^mir-Vli fi gcm'^ , (B5) 
aSSD = 1.7x 105(aSS^)-tmim§f"§/i gcm'^ , (B6) 



where 



/ = [i-(ri5££)^] 

r 



and risco is the radius of the innermost stable circular orbit. Note the dimcnsionless coordinate, r, is related to r by 
f = 2r. The parameter a describes the viscosity in the a-prescription (Shakura fc Sunvacv..l973.) . The superscripts 
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"ADAF" and "SSD" denote the Advection-dominated accretion flow and the Shakura-Sunyaev disk analytical solutions, 
respectively. To elucidate the general feature, we adopt typical values of the viscosity parameter such that a^^^ = 0.01 
and a^DAF ^ 

SOLVING THE RELATIVISTIC BERNOULLI EQUATION FOR THE LOCATION OF THE ALFVEN POINT 

Once the field strength, field geometry, and the mass loading per field line at different m are specified, we can 
calculate all the coefficients in the cold BE, Equation (26), and determine the location of the Alfven point, Ra- 

Since Equation (26) is solved in the geometrized unit (in length unit of GM/c^), extra transformations of units from 
c.g.s. unit are necessary, namely, 

H(gcom.) ^ B(c.g.s) X MbhG^/^/c" , (CI) 

Ai?y(geom.) - mpn-| (c.g.s) x MbhG^/'^/c^ = p J (c.g.s) x MbhG^/'^/c^ . (C2) 

These unit transformations can be done by two steps. First, we convert the unit from the c.g.s into the geometrized 
unit (Table 3). Second, we rescale the length unit in geometrized unit from "cm" into "GM/c^" (Table 4). 
Next, we find the location of the Alfven point Ra at various accretion rate by the following steps: 

1. Calculate the coefficients A and B in Equation (26). 

2. Guess Ra to compute E and L by Equations (29) and (30). 

3. Using E and L, calculate the coefiicients C and V in Equation (26). 

4. Using the calculated A , B, C, V, solve the cold BE equation(26). 

5. Find Ra so that the fluid can smoothly pass through the fast-magnetosonic point. 

A correct Ra will give a physical inflow solution that passes through the fast point before reaching the horizon 
with a nonzero speed (Figure 7). If Ra is incorrect, the resulting inflow solution becomes unphysical, being either 
sub- fast-magnetosonic (Figure 8), or discontinuous (Figure 9). 



TABLE 3 

Unit Conversion from c.g.s. into Geometrized Unit 



Quantity 


c.g.s Unit 


Geometrized Unit 


Dimension of the Geometrized Unit 


Rest-mass energy 


IJ, = rupc'' 


lx{c.g.s)G/c* = mp(c.g.s)G/c^ 


[cmj 


Number density 


n 


n(c.g.s) 


[cm-3] 


Magnetic field 


Bp 


Bp(c.g.s)Gi /c'^ 


[cm-l] 


Velocity 


Up 


up{c.g.s)/c 


[1] 


Mass loading 


7) = nUp/Bp 


r;(c.g.s)c/G2 


[cm-2] 



TABLE 4 

Converting the Length Unit from "cm" to "GM/c^" 


Quantity 


Geometrized Unit 


Geometrized Unit Scaled by GM/c** 


Dimension after Rescaled 


Rest-mass energy 
Number density 
Magnetic field 

Velocity 
Mass loading 


mp(c.g.s)G/c^ 
n(c.g.s) 

Bp{c.g.s)GVc2 
«p(c.g.s)/c 
»?(c.g.s)G/c2 


[mp{c.g.s)G/c^]/[(GM/c-)] 

[n{c.g.s)]/[{GM/c^)--^] 

[Bp(c.g.s)G5/c2]/[(GAf/c2)-i] 
up(c.g.s)/c 
[77(c.g.s)G/c2]/[(GM/c2)-2] 


[(GM/c-)] 

[(GM/c2)-3] 
[(GiU/c2)-i] 
[1] 

[(GM/c2)-2] 
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Fig. 7. Solution of the BE, Equation (26), when M/Mq = 10, e = 0.2 , rh = 0.001, and Ra/(GM/c^) = 6.128. The BE is a fourth-order 
algebraic equation and it has four roots at a fixed R in general. However, only the real solutions are shown here. This physical solution 
(I— >A— >F— >H) shows that the fluid is injected at the injection point (I) and passes through the Alfvcn (A) and the fast-magnetosonic (F) 
points, finally reaches the black hole (H) with a super-fast-magnetosonic speed. The vertical dotted line indicates the location of the inner 
light surface. We are not interested in the outward solution (— u'" < 0), because it represents an MHD flow escaping from the BH. 




Fig. 8. — Same as Figure 7, but for Rp^j (GM j (?) = 6.109. The solution is unphysical because the trajectory starting from the injection 
point reaches the horizon with a vanishing inward velocity. 
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R/(GM/c^) 



Fig. 9. — Same as Figure 7, but for Ra/{GM/c?) = 6.138. The solution that passes the Alfven point is unphysical, because the flow 
cannot reach the horizon. 



